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Abstract 

The dynamics of time-dependent coupled oscillator model for the 
charged particle motion subjected to a time-dependent external mag- 
netic field is investigated. We used canonical transformation approach 
for the classical treatment of the system, whereas unitary transformation 
approach is used when managing the system in the framework of quan- 
tum mechanics. For both approaches, the original system is transformed 
to a much more simple system that is the sum of two independent har- 
monic oscillators which have time-dependent frequencies. We therefore 
easily identified the wave functions in the transformed system with the 
help of invariant operator of the system. The full wave functions in the 
original system is derived from the inverse unitary transformation of the 
wave functions associated to the transformed system. 
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1 Introduction 



The time-dependent harmonic oscillators have attracted considerable in- 
terest in the literature thanks to their usefulness in describing the dy- 
namics of many physical systems. After the Bateman'spQ proposition 
concerning the use of time-dependent harmonic oscillator model in de- 
scribing dissipative systems, much attention was paid to quantum be- 
havior of nonconservative and nonlinear systems. 

In the meantime, coupled oscillators have emerged to become pow- 
erful modeling tools and, consequently, are frequently used in modeling 
wide range of physical phenomena. With the progress of research, one 
may be interested in what would happen if two-dimensional harmonic 
oscillator is elaborated through the coupling of two additive potentials? 
As far as we know, dealing with such an issue was set thirty years ago by 
Kim et al. [2-5]. Abdalla demonstrated how to treat the time-dependent 
coupled oscillators in the context of quantum mechanics[6]. The propa- 
gator for a time-dependent coupled and driven harmonic oscillators with 
time-varying frequencies and masses is investigated by Benamira (7j us- 
ing path integral methods. 

Among various systems that can be modeled by time- dependent cou- 
pled oscillators, the dynamics of charged particle motion in the presence 
of time-varying magnetic fields has played an important role in con- 
densed matter physics and plasma physics. There are plenty of applica- 
tions for this system such as magnetoresistance[5], the Aharonov-Bohm 
effect [9], magnetic confinement devices for fusion plasmas [TO], electro- 
magnetic lenses with variable magnetic fields [TT], cyclotron resonance |1 2], 
and entanglement of a two-qubit Heisenberg XY model[13]. Though all 
of these problems are interesting, we can find their exact analytic so- 
lutions only for a few special cases due to their complex mathematical 
structures. 

The quantum properties of a free electron, which have a time- 
dependent effective mass under the influence of external magnetic field, 
are investigated in both the Landau and the symmetric gauges [TU ITB] . 
Laroze and Rivera [TB] studied the dynamical behavior of electrons in the 
presence of a uniform time-dependent magnetic field and they presented 
the time evolution of the corresponding wave functions for the case that 
the initial state is a superposition of Landau levels. The propagators 
of a charged particle subjected to a time- dependent magnetic field are 
studied using the linear and the quadratic invariants |17j . 

Kim et al. [2-5] proposed a problem that what actually would take 
place if two harmonic oscillators are coupled so that the potential be- 
comes V(Xi,X 2 ) = \ {c.\X\ + C2X2 + C3X1X2) where C3 is a coupling 
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constant. They studied the corresponding density matrix in order to 
establish the Wigner function. In this work, we are interested in the 
problem of Hamiltonian that involves the coupling term X\X 2 in the 
presence of magnetic field. This system can be regarded as the general- 
ization of the Hamiltonian model given in Refs. [TT] and [18] . Though 
the coupling of two or more oscillators is among the most basic concepts 
in dealing with gyroscopic motions, interactions, and complex struc- 
tures, the related theory has been scarcely developed so far. This class 
of coupled harmonic oscillators can be used to describe numerous phys- 
ical systems. Some of them are the Bogoliubov transformation model of 
superconductivity [19], two-mode squeezed light [20], and the Lee model 
in quantum field theory [2T]. One of the main focuses of research carried 
out by Zhang et al. in connection with time-dependent coupled oscilla- 
tors including X\X 2 term are some specific problems of time-dependent 
coupled electronic circuits [22] 123]. 

We will use the invariant methods[24, 25] in order to derive the ex- 
act wave functions for time-dependent coupled oscillators in a variable 
magnetic field. The invariant operator method in describing the quan- 
tum features of time- dependent harmonic oscillators is firstly introduced 
by Lewis [21] and now became a very useful tool in developing quantum 
theory for the case where the Hamiltonian of the system is explicitly 
dependent on time. 

In Sec. 2, we formulate our problem by introducing a general time- 
dependent Hamiltonian describing the complicated motion of a charged 
particle in the presence of an arbitrary time-dependent magnetic field. 
Classical treatment of the system is presented in Sec. 3, on the basis of 
the canonical transformation method. Quantum analysis of the system 
is carried out in Sec. 4 using unitary transformation approach. The 
unitary transformation enables us to transform the original Hamiltonian 
(that is somewhat complicated) to that of a more simple system such as 
ordinary harmonic oscillator. We derive the quantum solutions of the 
system in Sec. 5 starting from the invariant operator associated to the 
transformed system described in Sec. 4. Finally, we give concluding 
remarks in the last section. 

2 Formulation of the problem 

For the dynamical system of our interest, the Hamiltonian has the form: 

H( Xl , X 2 , t) = ^y+^y+i m)Xl + C 2 (t)Xl + C^X^) , 

(1) 



3 



where Hi and n 2 are the conjugate momenta. Note that Hi and n 2 
can be simplified by choosing an appropriate gauge. Actually, in the 
symmetric gauge with ll( ~ B 2 ^ X 2 , ^p-X\, 0), they are given by 

n^-ffU.iW. + ffW (2) 

The parameters m x (t), m 2 (t), Ci(t), C 2 (t), and Cs(t) are arbitrary func- 
tions of time, (X X ,X 2 ) are the pair of position variables, and (Pi,P 2 ) 
are the canonical conjugate momentum variables. 

The main difference of our study from that of Ref. [16] is that we 
considered the coupling term XiX 2 in the Hamiltonian. Regarding the 
expressions of Hi and Il 2 , the Hamiltonian in Eq. (1) can be recasted 
into 



H(Xi, X^^ + ^ + l Ht)Xl + c 2 it)Xl + c 3 (t)XiX 2 ) 
+i (u 2c {t)P 2 X 1 - u Xc {t)P x X 2 ) , (3) 
where the new time-dependent functions c x (t), c 2 (t) and c 3 (t) are read 

ei(t) = C 1 (t)+m a (t)'£&, c 2 (t) = C 2 (t)+mi(t)^-, c a (t) = C 3 (t), 
with the cyclotron frequencies 



(4) 



, . eB(t) , , eB(t) 

mi(t) m 2 [t) 

3 Classical treatment 

The time-dependent canonical transformation approach is in fact very 
powerful in investigating the properties of dynamical systems described 
by a time-dependent Hamiltonian. In many cases, we can convert a 
given Hamiltonian into a simple and desired one by means of the canon- 
ical transformation. Therefore, in order to recast the solutions of this 
problem into a more soluble form, it is convenient to use the canoni- 
cal transformation method. To simplify the Hamiltonian given in Eq. 
(3), let us transform the variables (Xi, X 2 ,Pi, P 2 ) to the new variables 
(xi,x 2 ,p xl ,p X2 ) such that 

\m 2 (t)J \mx{t)J 
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Wi(w v«w)/ 

Replacing all of the canonical variables in Eq. (3) with the above ones, 
we have 

1 i 

H{x x , x 2 , t) = — — (p 2 Xi + p 2 X2 ) + - (di(t)xi + d 2 (t)xl + da(t)xix 2 ) 

+ ^fW 2 -^,), (8) 
where d\ — d% are new time-dependent functions of the form 

*w = C2(f) V2 = ( C2(i) + i-w^w) (=§) 1/2 • < 10 > 

<W«) = cs(t) = G 3 (t), (11) 

1 /2 

with the unique mass m(t) = (mi(t)m 2 (t)) and the cyclotron fre- 

1 /2 

quency u c (t) = (u lc (t)u 2c (t)) = eB(t)/m(t). 

To simplify the Hamiltonian of Eq. (8), we perform the following 
canonical transformation 



xA _/ cos0(t) sin0(t)\ (q\ 

PsA = ( cos0(t) sin0(t)\ /pA 
VxJ ~ V-sin0(t)cos0(t)y vp 2 y ' 



(12) 
(13) 



where 



<f>{t) = -~ / w c (*)d*. (14) 



If (gi, §2, V11V2) are canonical coordinates, there should exist a new Hamil- 
tonian H(qi, q 2 , t) which is determined by only in terms of the Hamilto- 
nian given in Eq. (8) with the aid of the linear transformation shown in 
Eqs. (12) and (13). The variables (xi,x 2 ,p xl ,p X2 ) and q 2 -,Px-,p 2 ) in 
two representations must satisfy the following relation 



(Pi<7i + P2<?2 - H(q 1 , q 2 , t) = p Xl Xi + p X2 x 2 - H(x 1 ,x 2 , t) + — — , (15) 

at 

where F\ is a time-dependent generating function in phase space, which 
should be determined afterwards. 
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From the fundamental equations known in classical mechanics 

d d 
Px 1 = Q^Fi(x 1 ,x 2i pi,P2,t), q 1 = —F 1 (xi,x 2 ,pi,p 2 ,t), (16) 

d d 
Px 2 = -Q^F 1 {x 1 ,x 2 ,pi,P2,t) , q 2 = —F 1 (x 1 ,x 2 ,pi,p2,t), (17) 

the generating function associated with the transformation is found to 
be 

F 1 (xi,x 2 ,p 1 ,p 2 t) = (pi cos0 + p 2 sin0) x\ + {—pi sin + p 2 cos 0) x 2 , 

(18) 

dFi ■ vj c (t) 

— = -(pit) {xxp X2 - x 2 p Xl ) = — [xip X2 - x 2 p xi ) . (19) 

In terms of the new conjugate variables (gi, q 2 ,Pi,p 2 ), the Hamiltonian 
of Eq. (8) becomes 

H ^ q2 ^ = 2^ft) ( p * +p i) + \( Xl ^ q * + + ^(^2) > ( 2 °) 

where 

Ai(i) = d\{t) cos 2 + d 2 {t) sin 2 — d 3 (t) sin0cos0, (21) 

\ 2 (t) = d 2 (t) cos 2 + d±(t) sin 2 + d 3 (t) sin cos 0, (22) 

A 3 (t) = 2 (di (t) - rf 2 (t) ) sin cos + d 3 (t) (cos 2 - sin 2 0) . (23) 

To eliminate the coupling term qiq 2 , we now perform the following canon- 
ical transformation |7J [22j [23] 



qi\ 1 / cos^f sin^\ /Q 

J " v^(t) Usin^cos^Hg 



(24) 



PA cosf sinf \ (PA \ (qi 

where 9(t) is an arbitrary function of time. Note that Eqs. (24) and 
(25) do not always represent the canonical transformation [27] between 
variables (qi,Pi) [i = 1,2] and (Qi,Pi). If (Qi,Pi) are canonical coordi- 
nates, there should exist a new Hamiltonian which is determined only 
by the Hamiltonian of Eq. (20) and the linear transformation given in 
Eqs. (24) and (25). The relation between variables (qi,Pi) and (Qi,Pi) 
in the two representations are 



2 2 Qp 

PiQi ~Hq = ~H q +—, (26) 

i=i i=i 
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where F is an another time-dependent generating function in phase 
space. 

Using the basic equations 
Vl = q 2 , P u P 2 , t), Qi = -^F (q u q 2 , P u P 2 , t) , (27) 

where i = 1,2, we see that the generating function is given by 
F (<?i, q 2 , Pi, P 2 , t) = ^~t) [p x cos ^ + P 2 sin q x 

+ ^Mt) (-P x sin ^ + P 2 cos q 2 
~m(t) {q\ + qf) . (28) 

Then, in terms of the new conjugate variables (Qi, Pi), the Hamiltonian 
can be represented in the form 

h q (Qu Q*, t) = \ (A 2 + Pi) + \n\{t)Q\ + \u? 2 {t)Q\ 

+ 9 ~Y- [P1Q2 - P2Q1] + 8{t)QiQi- (29) 
Here, the time- dependent coefficients Qi(t),Q 2 (t) and 5(t) are given by 

sm = (^ (t) cos , f + m S11 , m _ 1/2 , ,30) 
m) = (, ?(t) sin * f + ^ (t) cos * m + m*m) 1/2 , (31) 

S(t) = \ (m - Ol(t)) sm»( t ) + A3(t ^ 9(t) , (32) 



_ 2 Ai(t) { 1 (m 2 {t) m(t) 



where 



Qf(t) = - t rr + - ^J7T-2— tt , (33) 

1W m(t) 4 v m 2 (t) m(t)y' V ; 

^ w = w + i('^w_ 2 *(fly (34) 

2W m(f) 4 V^ 2 (t) m(t)y v ; 

If we take the choice 9{t) = Const, the terms P1Q2 and P2Q1 in Eq. 
(29) are canceled out so that the Hamiltonian becomes 

H Q (Q u Q 2 ,t) = l - (p? + pi) + ±ni(t)Qi+ 1 -nj(t)Qj + 5(t)Q 1 Q 2 . (35) 



Notice that, with the above canonical transformation, the coupling 
S(t) is a functional on the parameters of the original system. It is hence 
clear that the separation of variables in Eq. (35) requires that S(t) = 0, 
i.e. 

A 3 (t) = (Cul{t) -w?(t))m(*)tan0, (36) 

and consequently 

tang = - - ._^ 3 ^ _ 9 . .. . (37) 

By taking into account Eq. (36), the Hamiltonian in Eq. (35) is 
rewritten as 

H Q (Q 1 ,Q 2 ,t) = l - (Pi + Pi) + ]fli{t)Ql + ^n 2 2 (t)Ql (38) 

Then, Eq.(38) represents the sum of two independent Hamiltonians of 
the simple harmonic oscillators with the time- dependent frequencies 
and fl 2 (t). 

4 Quantum treatment 

The canonical transformations in classical mechanics, treated in the pre- 
vious section, is the analogous of the unitary transformations in quan- 
tum mechanics. Now we are going to demonstrate this relationship be- 
tween the two transformations and confirm how to obtain the quantum- 
mechanical Hamiltonian from the classical one. To manage the system 
in the context of quantum physics, we replace the canonical variables 
(Xi,X 2 ) in Eq. (3) by quantum operators (Xi,X 2 ). Then the corre- 
sponding Hamiltonian has the form 



+i (u 2c (t)P 2 X 1 - Wio(t)PiX 2 ) . (39) 

In this quantum case, the pair of momentum operators are given by 
(Pi = —ihd/dXi, P 2 = —ihd/dX 2 ). The Schrodinger equation in the 
original system is 

th^(X 1 ,X 2 ,t) = H(X 1 ,X 2 ,t)^(X 1 ,X 2 ,t). (40) 

To simplify the Hamiltonian in Eq. (39), we perform the unitary 
transformation such that 
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*(X 1 ,X 2 ,t) = U 1 (t)i/>(X u X 2 ,t), 
where U\{t) is a time-dependent unitary operator of the form 



(41) 



U\{t) =exp 



2h 



x exp 



2h 



+ X1P1) In 
(P 2 X 2 + X 2 P 2 ) In 



1/4' 



rai(t) 
m 2 (t) 

m 2 (t) N 1/4 
mi(i) 



(42) 



In this case, the Hamiltonian, Eq. (39), can be rewritten as 

1 



Hi(Xi,X 2 ,t) = 



2m(t) 



P1+P2 



+ 



+ 



X - [d x {t)X\ + d 2 (t)X\ + d z {t)X x X^ 



u c (t) 



P 2 X\ — P\X 2 



(43) 



It is easy to confirm that the commutation relations, [L z , X\ + 
Xf ] = and [L z , Pp + P 2 2 ] = 0, are hold where L z is the angular mo- 
mentum operator. This implies that there are common eigenfunctions 
between Lz and Xf + Xf , and between Lz and Pp + P|. However, 
Lz does not commutes with X\X 2 - [L^ , X\X 2 \ 7^ 0, and consequently 
[L z , P] 7^ 0. If we regard that L z and P do not have the same eigen- 
functions, it is not possible to simplify the Schrodinger equation 



at 



(44) 



by decomposing it. However, we can overcome this difficulty through 
the transformation of the Hamiltonian of Eq. (39) into a simple form 
by introducing an appropriate unitary transformation operators. In the 
first step, we perform the following unitary transformation 



il>(X 1 ,X 2 ,t) = U 2 (t)<p(X 1 ,X 2 ,t), 



(45) 



where 



P 2 (t) = exp(-— (P.Xx-PxX, 



w c {t)dt 



exp 



%L ; 



2h 



w c {t)dt 



(46) 
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Under this transformation, the Schrodinger equation (41) is mapped into 



th^(X u X 2 ,t) = H 2 (X 1 ,X 2 ,t)i P (X 1 ,X 2 ,t), 
at 

where the new Hamiltonian H 2 (Xi, X 2 ,t) has the form 

1 



(47) 



H 2 (Xi, X 2 ,t) 



2m(t) 



Pi + Pi 



+ - ( \ x {t)X{ + \ 2 {t)X l 2 + A 3 (t)XxX 2 



)■ 



(48) 



Now the term involving Lz has disappeared in Eq. (48). This means 
that the magnetic field term is removed in the new frame rotating with 
the time-dependent phase 4>{t) — — \ J zu c (t)dt. 

To decouple the Hamiltonian of Eq. (48), we take another unitary 
transformation such that 



<p(X u X 2 ,t) = V(t)x(X u X 2 ,t), 
where the unitary operator V(t) is given by 

V(t) = Vi{t)V 2 {t)%{t), 



(49) 
(50) 



with 



V 1 (t) = exp—[(P 1 X 1 +X 1 P 1 



x exp — 



)ln v / m(0] 
(P 2 X 2 + X 2 P 2 ) In ^Mt) 



V 2 (t) — exp 



-^(P^-Ax.) 



V 3 (t)=exp-— rh(t)(X* + Xt 



(51) 
(52) 

(53) 



Some algebra with the substitution of Eqs. (48) and (49) into Eq.(47) 
yields a transformed Hamiltonian that represents the sum of two uncou- 
pled simple harmonic oscillators having frequencies Qi(t) and Q 2 (t) and 
the unit mass: 



^ 3 (X 1 ,X 2 ,t) = V- 1 (t)i7 2 (X 1 ,X 2 ,t)V(t)-^- 1 (t)^(t) 

= \ (pi + H) + \n\{t)xl + \^ 2 {t)xl 



(54) 



At this stage, it is possible to confirm that the classically transformed 
Hamiltonian given in Eq. (38) is right, since the above equation is 
consistent with it. Note that U\{t) and Vi(t) given in Eqs. (42) and (51) 
are the squeeze operators, whereas U 2 (t) and V 2 (t) given in Eqs. (46) 
and (52) are the rotation operators characterized by the time-varying 
angles <p{t) and ^p, respectively. 
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5 Quantum solutions 



It can be seen that there exists invariant for the harmonic oscillator with 
time-dependent mass and/or frequency |24j. In our case, the transformed 
system consists of the two independent harmonic oscillators which have 
time-dependent frequency. It is easy to verify, from Liouville-von Neu- 
mann equation for the invariant I 



di di 1 rf . . 
-dl = m + Jh [I > H3] 



o. 



(55) 



that the invariant associated to the Hamiltonian of two-dimensional har- 
monic oscillator is given by 



I{X 1 ,X 2 ,t) = I{X 1 ,t) + I{X 2 ,t) 

2 



1 

2 



1 

+ 2 



) + ( Pl X x - p x X x 



I + I P2^2 - f>2 X 2 



(56) 



where pi(t) and P2{t) are c-number quantities obeying the auxiliary equa- 
tions 



p\ + n 2 1 (t) Pl = i/pl (57) 

p 2 + Ql(t)p 2 = l/pl (58) 

To guarantee the Hermiticity of Eq. (56) (F — I), we choose only the 
real solutions of the above two equations. It is clear that I(Xi,X 2 ,t) 
satisfies the Liouville-Von Neumann equation. We now derive a complete 
orthonormal set of eigenfunctions £ mn2 (Xi,X 2 , t) of I(Xi, X 2 , t) form the 
eigenvalue equation 

I(Xi, X 2 , t)£ nin2 (Xi, X 2 , t) = \ nin2 ^ niri2 (Xi, X 2 ,t), (59) 

where A„ in2 are time- independent eigenvalues. Through a straightfor- 
ward evaluation after inserting Eq. (56) into the above equation, we get 
the eigenvalues and the eigenfunctions such that 

A„ iri2 =h(ni + -j+h (n 2 + - J , (60) 
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Using Eqs. (42), (46), (50) and (65), we derive the full wave functions 
in the form 

1/2 



(X 1: X 2 ,t) 



^m x m 2 



x H, 



7r^ni!n 2 !2 n i+ n2 piP2j 

oi \ cos (4> + e/2)x 1 - 



W^sin (0 + #/2)X 2 



x H, 



n-2 



'mi sin 



h 1 / 2 Pl 

+ 0/2)X 1 + ym^cos 



+ 0/2) X 2 



x exp 



x exp 



imi 
27T 
irri2 
2h 



2 2 

2 + ^ 
2 2 



sin (0 + 20) 
sin (0 + 20) 



X\ 



xl 



x exp —Ajm x m 2 



7 ) cos(# + 20))X 1 X 2 



x exp« 



n\ + - 



dt' 

W) 



1 

n 2 + - 



dt' 



where the time-dependent coefficients j(t) and j3(t) are given as 

' pi i Id 
.Pi 



Pi 



2dt 



0{t) 



P2 



1 d 



p 2 P 2 2 ^ 



(68) 



(69) 



(70) 



The full solutions in the original system, given in Eq. (68), are exact 
since we did not use approximation or perturbation methods. Though 
these solutions are somewhat complicated, they are very useful in pre- 
dicting the quantum behavior of the system. A merit of such analytical 
solutions is that they can be employed in deriving the evolution of the 
probability distribution, regardless of the change of the system's param- 
eters. However, the numerical solutions in this field, such as the one 
obtained from FDTD (finite difference time domain) method [28]. are 
somewhat inconvenient as inputs to further analyses, since one should 
recalculate the results whenever the parameters of the system changes. 
Using Eq. (68), one can easily take a complete description of the charged 
particle motion even when the parameters of the system vary from time 
to time provided that the classical solutions of Eqs. (57) and (58) are 
known. 

6 Conclusion 

We investigated the quantal problem of the time-dependent coupled os- 
cillator model associated to the charged particle motion in the presence 
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of time-dependent magnetic field. Though the behavior of charged par- 
ticle in magnetic field drew great concern in both quantum and classical 
view point, researches in this line are rather concentrated on static prob- 
lems that can be modeled by time- independent harmonic oscillator. 

The system we treated in this work is however a more generalized 
one. It is summarized as follows: 

(i) We supposed that the effective mass of the charged particle varies 
explicitly with time under the influence of the time-dependent magnetic 
field. If electrons or holes in the condensed matter interact with environ- 
ment or various excitations such as pressure, energy, temperature, and 
stress, their effective mass may naturally vary with time|14j. Moreover, 
the random changes of the external field in the heterojunctions and solid 
solutions give rise to the variation of effective mass in accordance with 
the fluctuation of the composition in the system [29] . 

(ii) We let the external magnetic field B(t) be an arbitrary function of 
time. Therefore, the application of our theory is not confined in a special 
system that has a specific class of time-dependence for B(t). In fact, we 
can apply it in wide range of practical systems with the flexible choice 
of the type of B(t). 

(iii) Our system is further generalized by adding a coupling term X\X 2 
in the Hamiltonian. 

Through these generalization, the system became a somewhat com- 
plicated one that is described in terms of time-dependent Hamiltonian. 
Since the treatment of the original Hamiltonian system is not an easy 
task in this case, we transformed our system to that of a much more 
simplified one using two different techniques. In the first one, we carried 
out canonical transformations in order to simplify the problem relevant 
to the original classical Hamiltonian given in Eq. (1). After the trans- 
formation, the Hamiltonian reduced to a simple form associated to two 
uncoupled harmonic oscillators that each have time-dependent frequen- 
cies fli(t) and Q 2 (t). In the second technique we used an alternative 
approach on the basis of the unitary transformation method. With the 
choice of unitary operators Ui(t), U<z{t) and V(t), the quantum Hamil- 
tonian (39) has been transformed to an equally simple one as that of the 
canonical transformation previously performed, but within the realm of 
quantum mechanics. 

Since the Hamiltonian in the transformed system is very simple, 
we easily constructed dynamical invariant operator I(Xi, X 2 ,t) asso- 
ciated to the transformed system, as given in Eq. (55). The eigenstates 
£ ni „ 2 (Xi, X 2 , t) of this invariant operator are represented in terms of 
the Hermite polynomial. The Schrodinger solutions x nin2 (Xi, X 2 , t) in 
the transformed system are the same as £ >mn2 (X 1 , X 2 ,t) except for the 
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time-dependent phase factor e ia "i™2W. From the inverse transformation 
of x nin2 pTi, X 2 , t) with the unitary operators, we derived the full wave 
functions (quantum solutions) in the original system [see Eq. (68)]. The 
quantum solutions are expressed in terms of p\ and P2 that are the two 
independent solutions of the classical equation of motion given in Eqs. 
(56) and (57), respectively. Even if we represented the quantum solu- 
tions in terms of the classical solutions associated with the transformed 
system, it is also possible to represent them in terms of the classical 
solutions associated with original system. The wave functions given in 
Eq. (68) can be used to investigate various quantum properties of the 
system such as the fluctuations of canonical variables, the evolution of 
quantum energy, and probability densities, even when the parameters of 
the system vary from time to time. This is the advantage of such an- 
alytical solutions over numerical solutions obtained, for example, using 
the FDTD method [28]. 
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